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1.  Scientific  Objectives 

The  initial  aim  of  the  programme  was  to 
demonstrate  three  qubits  using  Na  donors  in  Si 
MOSFETs  by  the  end  of  the  programme.  The  initial 
objectives  can  be  summarised  as  follows:- 

2. 1  Fabricate  a  single  electron  transistor  using  AFM 
lithography  on  the  gate  of  a  MOSFET.  This  is  the 
detector  for  the  qubits. 

2.2  Fabricate  a  MOSFET  with  a  determined  number 
of  Na+  donors  in  the  oxide. 

2.3  Fabricate  and  detect  the  positioning  of  a  single 
Na+  ion  in  a  predetermined  position  to  allow 
fabrication  of  qubits. 

2.4  Measure  the  T2  decoherence  time  for  Na  in  Si 
using  the  temperature  dependence  of  the  impurity 
band  and  spin  dependent  scattering. 

2.5  Provide  a  model  for  the  Na  donor  fabrication. 

2.  Summary  of  Important  Results 
2.1  Electronic  Transport 

Since  the  invention  of  the  silicon  MOSFET, 
understanding  the  influence  of  impurities,  espe¬ 
cially  sodium  contamination,  on  device  perform¬ 
ance  has  been  a  priority  and  continues  to  provide  a 
rich  system  for  investigation  by  experimental  and 
theoretical  physicists  alike.  The  electronic  proper¬ 
ties  of  sodium  doped  MOSFETs  were  first  studied 
by  Fowler  and  Hartstein  [1,2]  in  the  1970s.  They 


reported  a  single,  broad  peak  in  the  subthreshold 
drain  current  against  gate  voltage  and  attributed  it 
to  the  formation  of  an  impurity  band  induced  by  the 
presence  of  sodium  ions  near  the  Si-Si02  interface. 
Further  studies  of  narrow  channel  devices  (~100 
nm)  demonstrated  a  series  of  reproducible  sharp 
peaks  [3],  while  later  experiments  found  evidence 
for  resonant  tunneling  between  localised  states  in 
the  channel.  [4,5,6]  For  sufficiently  low  impurity 
concentrations,  the  overlap  between  neighbouring 
localized  electron  wavefunctions  and  consequently 
the  hybridisation  of  their  excited  states  is  predicted 
to  be  reduced  [7],  splitting  the  single  impurity  band 
observed  at  high  concentrations  into  the  ground  and 
excited  bands  as  modeled  by  Ghazali.  [8]  Increasing 
the  resistivity  of  the  silicon  substrate  reduces  the 
scattering  from  acceptors  at  the  Si-Si02  interface, 
allowing  the  possibility  for  such  a  band  splitting  to 
be  experimentally  observed  in  the  transport.  In  this 
paper,  we  will  present  evidence  for  the  observation 
of  two  separate  impurity  bands  with  a  soft  gap,  based 
on  analysis  of  the  temperature  dependent  conductiv¬ 
ity  below  20  K. 

The  device  we  used  is  a  MOSFET  fabricated 
on  a  (100)  oriented  p-silicon  wafer  and  was  subse¬ 
quently  patterned  in  the  circular  Corbino  geometry 
to  eliminate  Hall  voltages  and  possible  leakage 
paths.  The  effective  gate  channel  length  and  interior 
width  were  respectively  1  //m  and  346  //m.  A  high 
resistivity  wafer  ( 104  Q.cm)  provided  a  background 
concentration  of  less  than  10 12  cm  3  of  boron  corre- 
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sponding  to  a  mean  distance  between  impurities  of 
1  pim.  A  35  nm  gate  oxide  was  grown  at  950°C  in  a 
dry,  chlorine-free  oxygen  atmosphere.  The  phospho¬ 
rous  implanted  and  aluminium  sputtered  contacts 
were  highly  metallic  and  Ohmic  at  all  temperatures 
investigated.  Sodium  ions  were  introduced  onto  the 
oxide  surface  by  immersing  the  device  in  a  10'7  N  ? 
solution  of  high  purity  sodium  chloride  (99.999%) 
in  de-ionised  water. 

The  surface  of  the  chip  was  dried  with  nitrogen 
gas  and  an  aluminium  gate  subsequently  evaporated. 

To  observe  or  remove  the  low  temperature  conduc¬ 
tivity  structures,  the  mobile  ions  are  drifted  through 
the  oxide  to  the  Si/Si02  interface,  or  returned  to  the 
Al-Si02  interface  by  applying  either  a  +4  V  or  a  -4 
V  DC  gate- substrate  bias  for  10  min  at  65°C  before 
the  device  is  cooled  down  to  helium  temperature 
at  which  sodium  looses  its  diffusivity  in  the  oxide. 
All  measurements  were  performed  using  standard 
low-noise  lockin  techniques  with  an  amplifier  of 
108  V  /  A.  The  AC  excitation  was  maintained  at  15 
pi V  with  a  frequency  of  1 1  Hz.  Suitable  RC  filters 
were  employed  to  eliminate  any  DC  offset  from 
the  amplifier.  The  gate  voltage  was  controlled  by 
a  high  resolution  digital  to  analog  converter.  All 
experiments  were  performed  in  a  3He  cryostat  and 
the  temperature  was  measured  by  a  calibrated  ger¬ 
manium  thermometer. 

Fig.  1  shows  the  conductivity  a  of  our  device 
at  300  mK  versus  gate  voltage  Vg  for  the  case  where 
the  sodium  ions  had  been  drifted  to  the  Si-Si02  inter¬ 
face.  Two  groups  of  peaks  appear  clustered  around 
Vg  =  -2V  and  -0.5  V  and  separated  by  a  region  of 
low  conductivity  and  limited  by  noise.  The  origin 
of  the  peaks  themselves  will  be  discussed  later. 
Following  a  -4  V  drift,  no  structure  was  detectable 
over  the  full  range  of  gate  voltages  investigated  but 
a  difference  in  threshold  voltage  of  0.2  V  was  found 
at  77  K  between  the  characteristics  of  the  device 
following  +  4  V  and  -4  V  drifts.  This  was  attributed 
to  the  presence  of  mobile  charges  close  to  the  Si- 
Si02  interface  at  a  concentration  of  3.7x  1011  cm'2 
corresponding  to  a  mean  impurity  separation  of  16 
±  1  nm.  In  a  reference  device  where  no  sodium  was 
introduced,  no  sub-threshold  conductivity  peaks  and 
no  shift  of  the  threshold  voltage  appeared  for  any 
drift  conditions  investigated.  Tunneling  through  the 


Fig.  1:  The  source-drain  conductivity  versus  gate  voltage  at 
300  mK  following  a  +4  V  drift  and  a -4  V  drift  ( dotted 
line). 

oxide  and  other  leakage  currents  were  discounted  as 
the  gate  leakage  current  was  below  50  fA  at  4.2  K 
and  was  approximately  constant  over  the  range  of 
gate  voltages  used.  We  remark  also  that  no  hysteresis 
and  thus  no  charging  effects  that  have  been  reported 
in  similar  devices  [9]  were  observed  here. 

The  presence  of  two  distinct  ranges  of  Vg  where 
peaks  appear  suggests  the  possibility  of  a  split  impu¬ 
rity  band  (Fig.  1).  Such  a  splitting  into  a  ground  and 
an  excited  state  is  expected  to  happen  for  low-dop¬ 
ing  concentration  [8,10]  if  one  takes  into  account  the 
overlaps  between  impurity  wavefunctions  and  uses 
a  multi-band  formalism  [11].  The  conductivity  of  a 
Si-MOSFET  is  not  directly  related  to  the  density  of 
states  and  the  fact  that  we  see  two  regions  of  high 
conductivity  separated  by  a  region  of  low  conduc¬ 
tivity  is  only  an  indirect  indication  that  the  device 
density  of  states  consists  of  two  bands  separated 
by  a  gap.  In  fact,  through  the  Kubo  formalism,  con¬ 
ductivity  tends  to  be  related  to  local  paths  through 
a  disordered  device  but  density  of  states  is  a  global 
property.  In  order  to  show  that  the  density  of  states 
splits  into  two  bands,  we  have  looked  at  the  tempera¬ 
ture  dependence  of  the  conductivity  at  a  series  of 
different  gate  voltages.  For  all  gate  voltages  studied, 
the  conductivity  decreases  non-monotonically  as 
temperature  is  lowered  (Fig.  2).  In  the  range  1  K 
to  20  K,  we  observe  the  characteristics  of  hopping 
conduction  so  that  the  conductivity  g(T)  is  fitted  to 


Page  3 


DAAD19-01 -1-0552 


University  of  Cambridge 


*  <  ; 

!M»», 

3 ' 

a  • 

O  £  o 

o 

•  o 

<,  !.  5}/ 

vt** 

b) 

io-2 


0,4 

s 


0,5 


Fig .  2 


Temperature  dependence  of  the  conductivity  for  V  = 
0.1  (o),  -0.24  (•),  -1.48  (*)  and -2.26 (<\),  b)  variation 
of  the  reduced  %2  with  s  =  0.412,  0.406,  0.394  and  0.385 
respectively  and  c)  variation  of  the  reduced  %2  with  y 
for  s  equal  to  the  optimum  value  found  in  2(b),  for  the 
appropriate  gate  voltages,  displaying  consistency  with 
a  minimum  aty  =  2  ( dotted  line).  Lines  in  b)  represent 
the  exponent  for  the  Mott  hopping  regime  (left)  and 
Efros  Schlovskii  regime  (right). 

the  generalised  equation: 


<7 


=  a0T^e-m 


(1) 


where  T0  and  G0  depend  on  gate  voltage. 

The  best  values  for  y  and  s  are  found  by  mini¬ 
mizing  the  value  of  the  reduced  chi-square  deviation 
%2  with  standard  procedures.  Figs  2(  a)-  (b)  show  the 
optimum  values  determined  for  four  gate  voltages, 
one  point  in  the  threshold  region  (V  =0.1  V),  the 
upper  band  (Vg  =  -0.24  V),  the  gap  (Vg  =  -1.48  V) 
and  the  lower  band  (V  =  -2.26  V).  The  resulting 
fits  for  g(T)  are  valid  over  three  orders  of  magni¬ 
tude  in  G  (Fig.  3).  Studies  in  gate  voltage  show  that 
y  =  1.98  ±  0.04  and  s  =  0.39  ±  0.02  for  V  below 
0.25  V.  Above  this  point,  the  value  of  s  decreases 
rapidly  towards  s  =  1/3  for  T  >  4K  and  the  range  is 
better  described  by  Mott  hopping  conduction  [12]. 
It  is  worth  noticing  that  the  smallest  values  of  s  are 
obtained  for  band  centre  regions  -0.8  V  <  Vg  <  -0.5 
V  and  -2.2  V  <  Vg  <  -2.1  V  for  which  the  hopping 
lengths  are  relatively  smaller  and  the  Coulomb 
interactions  weaker. 

We  emphasise  that  the  use  of  temperature 
dependent  exponential  prefactors  allows  the  fine 
distinction  between  Mott  (s  =  1/3),  Efros-Schk- 
lovskii  (s  =  1/2)  [13]  and  the  regime  under  study 


Fig.  3:  Temperature  dependence  of  the  conductivity  for  V 
=  0.1  (of  -0.24  (•),  -1.48  (*)  and  -2.26  (< )  for  the 
optimum  values  of  g  and  s  defined  in  Eq.  1 .  For  clarity, 
curves  are  shifted  downwards  respectively  by  0,  1,  2 
and  3  from  their  original  values. 

(s  ~  0.39).  The  finding  that  y  ~  2  is  consistent  with 
the  formulation  for  g(T)  given  by  Allen  and  Adkins 
[14,15]  if  it  is  rederived  for  the  2D  case.  This  gives 
the  conductivity  in  terms  of  both  the  localization 
length  ^  and  the  density  of  states  at  the  Fermi  level 
n(Ep).  Accounting  for  the  fact  that  Coulomb  interac¬ 
tions  between  electrons  in  different  localized  states 
modifies  the  density  of  states  close  to  the  Fermi  level 
so  that  n(E)  =  N0 1  E  -  Ep  lp,  we  obtain: 


m  —  -i(  P+1  )  _(Zb.)p+3 

cr  =  cr0T  7U+3)e  V  t  )  with  y  =  2 

(2) 

<70  =  A0  \n2N^3p+13  (p+  lf3l/raC“1P‘ 


ksTo  = 


p  +  3  P+1 

p  +  1 


(p  +  l)3 

irNo^ 


1 

p+1 


(3) 

(4) 


where  A0  is  a  constant  depending  on  the  electronic 
properties  of  bulk  silicon  and  kB  the  Boltzmann 
constant. 

The  conditions  given  by  Allen  [14]  for  the  use 
of  the  equations  Eqs.  2-4  are  satisfied  in  our  device. 
Recent  calculations  showed  that  sodium  ions  in  the 
oxide  may  trap  either  one  or  two  valence  electrons 
against  the  Si/Si02  interface  and  that  the  wavefunc- 
tions  of  the  localized  states  remain  hydrogen-like. 
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[16]  Also,  hopping  conduction  is  present  in  our 
device  so  that  the  only  mobile  electrons  are  found  in 
an  energy  band  corresponding  to  a  hopping  energy 
around  the  Fermi  level.  The  resulting  Coulomb  gap 
for  s  =  0.39  (i.e.  n(E)  ~  I E  -  Ep  I0  30)  is  much  sharper 
than  for  the  Efros  regime  where  n(E)  ~  I  E  -  Ep  I. 
This  behaviour  has  been  predicted  by  Blanter  and 
Raikh  [17]  while  studying  2D  systems  localized 
by  disorder.  They  have  shown  that  a  metallic  gate 
close  to  the  interface  provides  image  potentials 
that  modify  the  densityofstates  to  the  form  n(E)  = 
n(Ep)+N0  IE-Ep  1 1/3  at  T  =  0  K  and  thus  gives  an  ex¬ 
ponential  dependence  of  T-0  4  for  the  conductivity. 
This  behaviour  is  explained  by  the  fact  that  the  oxide 
thickness  plays  the  role  of  the  screening  length  and 
that  initial  and  final  states  become  electrostatically 
independent  at  low  temperature  when  the  hopping 
length  R  becomes  greater  than  twice  the  oxide 
thickness  d.  This  then  produces  a  crossover  from  the 
Efros  to  the  Mott  regime.  The  localization  length  in 
our  device  was  approximately  estimated  by  fitting 
the  conductivity  with  p  =  0,  using  the  Mott  formula 
for  kBTQ  and  taking  the  2D  value  for  the  density  of 
states.  For  Vg  =  -0.4  V,  this  gives  £,  ~  20  nm  and  a 
hopping  length  R  ~  63  nm  at  1  K.  Our  device  has  a 
gate  oxide  of  d  =  35  nm  and  the  Coulomb  interac¬ 
tions  may  be  screened  by  the  electrostatic  gate  as 
R  ~  2d.  This  demonstrates  the  device  may  well  be 
in  the  regime  described  by  Blanter  and  Raikh.  Fits 
of  the  conductivity  using  y  =  2  and  s  =  0.4  are  still 
valid  for  V  smaller  than  0.25  V  and  from  1  K  to  1 8 

g 

K  typically,  so  we  will  proceed  with  these  values. 
In  this  regime,  the  density  of  states  at  Ep  is  given 
by  : 


n  (Ep)  =  no  (Ep) 


2irno  (Ep)2  de2 

47re0er  (5) 


no  (Ep) 


f  No\  ( 47reoer\  3 
v  2tt  )  V  2 d2d2  ) 


(6) 


where  er  is  the  permittivity  of  silicon. 

These  equations  have  been  derived  at  T  =  0 
under  the  condition  that  n(Ep)d2V(d)  (<  0.095  in  our 
device)  is  small  where  V  is  the  impurity  Coulomb 
potential.  At  higher  temperature  but  for  T  «  V (d),  the 
gap  is  partly  filled.  The  Coulomb  interactions  then 


become  negligible  for  energies  below  V(2d)  and  the 
density  of  states  saturates  such  as  n0(Ep)  =  n0(Ep  + 
V(2d)).[18]  This  condition  decreases  the  value  of 
v  from  1  down  to  0.11.  Combining  (3)  and  (4),  the 
parameters  N0  and  £,  were  expressed  in  terms  of  T0 
and  G0,  values  easily  accessed  experimentally.  The 
density  of  states  was  extracted  using  equation  (5) 
and  (6).  Several  values  of  S  are  found  in  the  litera¬ 
ture  [19]  but  the  most  recent  value  was  found  by 
Martin  [20]  who  gives  S  =  16.14eV  implying  A0  = 
1.10  x  10'3  ekh  '.K-0  8. 

The  density  of  states  is  shown  in  Fig.  4.  The 
value  we  find  is  two  orders  of  magnitude  lower 
than  the  pure  2D  case  (‘”  1.6  x  1014  eV'.cnr2  with 
valley  degeneracies).  The  large  background  is  pre¬ 
dominantly  due  to  the  conduction  band  tail  which 
spreads  over  the  full  range  of  Vg  studied.  The  up¬ 
per  band  also  has  a  significant  tail.  The  presence 
of  density  of  states  tails  is  a  common  occurance  in 
disordered  systems  with  low  impurity  concentra¬ 
tion  and  localized  wavefunctions  [21,22,23]  but 
its  linear  shape  for  Vg  <  0  V  is  unusual.  It  has  been 
attributed  to  the  formation  of  regions  of  constant 
local  potential  energy  at  the  Si-Si02  interface  and 
containing  a  random  number  of  charges.  [24]  Two 
regions  of  higher  density  are  superimposed  on  the 
background  density  and  correspond  to  the  upper 
and  lower  groups  of  peaks.  This  confirms  that  the 
structure  observed  in  the  conductivity  is  due  to  the 
pres-ience  of  two  separate  bands.  By  numerically 


The  dashed  and  dotted  lines  represent  respectively  the 
background  density  due  to  the  conduction  band  tail 
and  the  upper  band  tail. 
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Fig.  5:  Variation  of  the  localisation  length  (•)  and  the 
trnasmission  coefficient  at  1  K(o)  with  gate  voltage. 

subtracting  the  background  density  and  integrating 
over  the  appropriate  gate  voltage  and  by  supposing 
a  linear  relation  between  the  gate  voltage  and  the 
surface  potential  energy,  we  estimate  the  upper  band 
contains  approximately  3  times  the  number  of  states 
as  the  lower  band.  This  value  is  an  upper  bound  as 
energetically  deeper  states  could  not  be  accessed 
experimentally. 

The  variation  of  the  localization  length  £,  (Fig. 
5)  follows  that  of  the  density  of  states,  showing  the 
two  bands  do  correspond  to  the  more  conductive 
regions.  The  value  of  £,  decreases  rapidly  when  ap¬ 
proaching  the  threshold  voltage.  This  is  expected 
as  the  region  between  Vg  =  0  V  and  Vg  =  0.4  V  is 
both  in  the  band  tails  of  the  conduction  band  and 
the  upper  band,  the  conduction  band  edge  being 
well  above  0.4  V.  Thus,  this  region  is  still  a  region 
of  strong  localization.  Nevertheless,  the  value  of 
is  expected  to  rise  once  the  Fermi  energy  crosses 
the  conduction  band  edge.  Finally,  taking  into  ac¬ 
count  fitting  errors  as  well  as  the  discrepancy  in  the 
value  of  the  deformation  potential,  we  estimate  the 
localization  length  within  5  %  and  the  density  of 
states  within  6  %.  From  the  same  derivation  which 
gave  (1),  (2)  and  (3),  the  transmission  coefficient 
is  obtained: 


In  (Te) 


2  (To\p+3 

P  +  3  \TJ 


(7) 


The  experimental  variation  of  the  transmission 
coefficient  in  gate  voltage  (Fig.  5)  shows  that  the 


conductivity  of  the  two  bands  mostly  comes  from 
the  higher  mobility  of  the  states  of  energies  within 
the  upper  and  lower  bands  and  less  from  the  increase 
in  the  density  of  states. 

In  conclusion,  we  have  observed  in  electrical 
transport  measurements,  an  unusual  hopping  regime 
with  an  exponent  0.4  which  results  from  the  screen¬ 
ing  of  the  Coulomb  interactions  by  the  metal  gate. 
We  have  shown  that,  consequently,  both  the  localiza¬ 
tion  length  and  the  density  of  states  can  be  extracted 
from  the  temperature  dependence  of  the  conductiv¬ 
ity.  This  analysis  has  given  strong  evidence  for  the 
existence  of  two  separate  bands  and  a  soft  gap  at 
low  temperature.  This  may  result  from  the  splitting 
of  the  impurity  band  into  a  lower  and  an  upper  band 
in  presence  of  Coulomb  interactions.  The  formation 
of  the  two  bands  results  itself  from  the  presence  of 
a  low  concentration  of  sodium  impurities  close  to 
the  Si-Si02  interface.  The  conditions  for  the  obser¬ 
vation  of  such  a  formation  may  be  the  creation  of 
deep  but  well  separated  impurity  potentials  at  the 
interface,  resulting  in  electron  localization.  Finally, 
the  electron  screening  may  be  sufficiently  weak  and 
the  disorder  not  too  strong  to  allow  a  Mott-Hubbard 
transition  to  take  place.  The  two  bands  could  then 
possibly  be  Hubbard  bands. 

2.2  Electron  Paramagnetic  Resonance 

A  commercial  Briiker  EMX  spectrometer  was 
used  for  the  measurements.  To  increase  sensitivity, 
a  modulated  field  AB  .  is  applied  at  frequency  to  . 
in  addition  to  the  dc  microwave  power.  The  Briiker 
was  used  with  a  modulated  frequency  of  100  kHz 
and  the  modulation  amplitude  was  set  to  be  suffi¬ 
ciently  small  to  prevent  additional  broadening.  The 
lineshape  is  assumed  to  relax  between  successive 
passes  of  the  microwave  resonance  condition.  In 
the  case  of  homogeneous  broadening,  the  linewidth 
arises  solely  due  to  Tt  and  T2  relaxation  processes  at 
the  EPR  centres.  A  continuous  flow  cryostat  allows 
the  sample  inside  the  microwave  cavity  to  be  cooled 
down  to  about  10  K  while  radiation  of  frequency 

9.2  GHz  is  applied  and  measured  with  a  HP53 15 1 A 
frequency  counter. 

A  batch  of  MOSFETs  with  210  nm  gate  oxides 
was  prepared  and  dipped  in  10-5  M  NaCl  solution 
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Fig.  6:  A  MOSFET  capacitor  dipped  in  10~5  M  so¬ 
lution  ofNaCl.  A  -12.9  V  threshold  shift  is  ob¬ 
served  with  +10  V corresponding  to  a  1.32xl012 
cm~2  sodium  density. 

(1.32xl012  cm-2  Na  from  CV  of  Fig.  6).  The  deriva¬ 
tive  of  some  of  the  spectra  are  shown  in  Fig.  7  about 
6  days  after  the  Na  was  put  into  the  oxide.  The  intial 
spectrum  shows  two  phosphorus  peaks  symmetrical 
around  g=1.9992  and  split  by  41.94  Gauss  due  to 
the  hyperfine  interaction  with  the  31P  nucleus  [25]. 
This  must  be  phosphorus  in  the  substrate  at  ~1013 
cm-3  concentrations  even  though  the  substrates  are 
p-type  to  5xl014  cm-3.  Around  g~2  there  is  a  peak 
close  to  the  electron  g-value.  After  applying  +10 

V  the  peak  disappears  and  then  reappears  after  -10 

V  bias.  A  further  positive  bias  results  in  the  peak 
disappearing. 

Again  this  peak  may  be  due  to  electrons  from 
the  sodium  impurity  band  but  also  it  may  be  due  to 
defects  at  the  Si/Si02  interface.  Fig.  8  demonstrates 
the  diffusion  of  Na  in  the  gate  oxide  of  a  MOSFET 
after  the  6  days  from  fabrication.  If  a  diffusion 
constant  of  2  x  10~16  cm2/s  is  used  then  a  significant 
amount  of  Na  can  diffuse  to  the  oxide  interface. 
While  data  exists  in  the  literature,  there  is  a  large 
spread  in  the  diffusivities  with  data  ranging  from  2.8 
x  10-15  m2/s  to  10-19  m2/s.  Clearly  there  is  sodium  at 
the  Si/Si02  interface  after  6  days  if  the  higher  dif- 
fusivity  is  used  and  therefore  it  is  possible  that  the 
initial  measurement  in  Fig.  7  could  be  from  a  sodium 
impurity  band.  The + 10  V  would  take  the  MOSFET 
into  strong  inversion  with  a  carrier  density  of  2.6 
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Fig.  7:  The  derivatives  of  4  spectrum  taken  on  a 
single  sample  for  the  last  run  at  10  K  (off  set 
vertically  for  clarity).  The  two  peaks  symmetri¬ 
cal  around  g= 1.9992  are  from  P  impurities.  The 
central  peak  is  close  to  the  free  electron  value 
of  g= 1.99869. 


x  1012  cm-2  which  is  over  twice  the  Na  impurity 
density  measured  by  CV.  While  EPR  resonances 
from  inversion  layers  have  been  observed,  they  are 
typically  a  sharp  downward  peak  in  the  derivative 
with  two  smaller  satellite  peaks  [25],  These  inver- 


Distance  from  oxide  surface  (nm) 


Fig.  8:  The  concentration  of  Na  in  the  gate  oxide 
of  the  MOSFET  as  a  function  of  distance  from 
the  Aluminium  gate  for  4  diffusivities.  Even  at 
the  Si/Si02  interface  which  is  210  nm  below  the 
surface,  a  finite  amount  of  sodium  has  diffused 
after  6  days. 
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sion  layer  peaks  are  clearly  not  observed  at  +10  V,  0 
V  or  -10  V.  The  threshold  voltage  for  this  device  is 
-10.4  V.  Therefore  measuring  at  -10  V  would  have 
the  Fermi  energy  just  above  the  sodium  impurity 
band  but  without  any  significant  formation  of  a  high 
density  inversion  layer. 

The  calculated  wavefunction  has  a  p-like  or¬ 
bital  and  therefore  there  can  be  no  hyperfine  inter¬ 
action  with  the  nuclear  spin  of  the  I=+3/2  sodium 
nucleus.  If  a  hyperfine  interaction  was  present,  four 
resonances  would  be  observed.  As  the  hyperfine 
interaction  is  the  limiting  factor  for  T2  in  phospho¬ 
rus  doped  Si,  this  may  be  beneficial  for  the  present 
sodium  system.  This  also  means  that  no  hyperfine 
splitting  of  the  any  peak  from  electrons  in  the  so¬ 
dium  impurity  band  should  be  expected  from  inter¬ 
actions  with  the  I=+3/2  nuclear  spin  of  23Na  which 
is  the  most  abundant  isotope  of  sodium. 

It  is  also  possible  that  the  resonance  in  Fig.  7  is 
due  to  defects  at  the  Si/Si02  interface  [26-28],  The 
main  interface  defect  peak  is  due  to  the  pb0  which 
has  a  g- value  of  2.006  for  the  [100]  orientation  of 
the  present  sample.  Fig.  9  shows  the  integrated  lines 
from  Fig.  7  and  the  pb0  peak  can  be  observed  at  2.006. 
The  g-values  for  this  graph  have  been  calibrated  to 
the  correct  positions  using  the  two  phosphorus  peaks 
which  are  symmetrical  around  g=  1.9992  ±  0.0494. 
The  pb0  peak  is  therefore  the  broad  resonance  around 
2.006  and  a  sharp  resonance  remains  at  1 .998  with 
a  T2*  of  142  ns.  The  only  other  known  defect  state 
is  the  E’  centre  [29]  which  is  centred  at  g=2.0011. 
This  peak  has  been  shown  to  saturate  at  /iW  pow- 


g  Value 

Fig.  9:  The  integrated  lines  from  Fig.  12  with  a 
number  of  curves  fitted  to  the  data. 


ers  and  the  present  measurements  demonstrate  the 
peak  at  mW  powers  where  the  defect  should  not 
be  observed.  In  addition  extreme  measures  such  as 
high  temperature  anneals  in  02  or  y  radiation  had 
to  be  used  to  form  such  defects.  All  this  evidence 
strongly  suggests  that  the  sharp  peak  around  1.998 
is  due  to  electron  spins  from  the  sodium  impurity 
band.  The  only  other  posibility  is  that  an  unknown 
Si/SiO?  impurity  state  due  to  charge  in  the  oxide 
might  create  some  line  at  this  particular  value.  No 
such  states  have  been  observed  before  and  if  such 
a  state  was  created  it  would  probably  be  the  result 
of  the  sodium  in  the  oxide. 

Fig.  10  shows  the  Vogt  function  fitted  to  the 
potential  sodium  impurity  band  peak  at  g=  1.998. 
From  the  width  of  this  fit,  the  inverse  gives  the  lower 
bound  for  the  T2  lifetime  [25]  to  be  0.5  /is  (see  Fig. 
11  for  a  description  of  T2  and  T  *).  Recently  new 
pulsed  EPR  measurements  of  phosphorous  doped 
silicon  have  appeared  in  the  literature  [].  Using 
isotopically  purified  silicon  with  a  phosphorus 
concentration  of  1015  cm-3  a  T2  lifetime  of  14  ms 
has  been  measured  at  7  K.  One  potential  issue  is 
that  the  present  technique  being  employed  cannot 
detect  T2  lifetimes  longer  than  100  /is.  Therefore  if 
T2  for  electrons  in  the  sodium  impurity  band  is  as 
long  as  most  of  the  other  shallow  donors  in  silicon, 
the  present  technique  cannot  measure  it  unless  T2 

- Response 

200  "I  . . Fit  of  Voigt  function 
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Fig.  10:  The  Voigt  function  fitted  to  the  potential 
sodium  impurity  band  peak  at  g= 1.998  to  ob¬ 
tain  Tr  The  extracted  T2from  this  is  >  0.5  psfor 
electrons  5.2  nm  apart. 
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Lorentzian  Gaussian 


Fig.  11:  The  difference  between  homogenously  and 
inhomo genously  broadened  EPR  peaks  and  the 
extraction  of  the  appropriate  lifetime. 

is  reduced  by  moving  to  higher  temperatures.  High 
temperature  measurements  were  also  taken  for  the 
last  set  of  samples,  the  phosphorus  peaks  were  ob¬ 
served  up  to  high  temperatures  but  no  other  peaks 
became  observable  at  higher  temperatures  around 
g=  1.99869. 

In  summary,  in  two  sets  of  measurements  we 
have  achieved  a  signal  which  is  close  to  the  expected 
resonance  at  g=  1 .99869  and  obtained  a  lower  bound 
for  the  T2  lifetime  from  these  resonances  of  0.5  ps 
for  a  sodium  density  of  1.32  x  1012  cm-2.  The  state 
is  not  hyperfinely  split  as  expected  from  a  p-orbital 
wavefunction  and  does  not  correspond  to  any  of  the 
known  interface  defect  resonances.  The  spacing  of 
electrons  for  the  lowest  density  we  have  managed 
to  measure  is  only  5.2  nm  which  is  still  close  to 
the  Bohr  radius  and  hence  exchange  interactions 
between  the  electrons  will  still  be  quite  strong  and 
therefore  reduce  the  T2  lifetime.  For  qubits,  this 
electron  spacing  will  be  significantly  larger  and 
the  exchange  interaction  should  be  much  smaller. 
Also  it  is  the  Orbach  process  which  typically  limits 
the  T2  lifetime  in  semiconductors  and  this  has  a  T7 
temperature  dependence  [30].  Therefore  operating 
a  quantum  computer  below  1  K  should  allow  sig¬ 
nificantly  higher  T2  lifetimes  to  be  achieved  than  the 
0.5  ps  measured  in  the  present  results. 
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